|L. N. Bulaevskii and S. Z. Lin, Phys. Rev. B 86. 224513 (2012).| 



Self-induced pinning of vortices in the presence of ac driving force in magnetic superconductors 
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We derive the response of the magnetic superconductors in the vortex state to the ac Lorentz force, 
Fi(t) = sin(ti;f), taking into account the interaction of vortices with the magnetic moments described by 
the relaxation dynamics (polaronic effect). At low amplitudes of the driving force Fac the dissipation in the 
system is suppressed due to the enhancement of the effective viscosity at low frequencies and due to formation 
of the magnetic pinning at high frequencies a>. In the adiabatic limit with low frequencies ai and high amplitude 
of the driving force F„^, the vortex and magnetic polarization form a vortex polaron when Fi(t) is small. When 
Ft increases, the vortex polaron accelerates and at a threshold driving force, the vortex polaron dissociates and 
the motion of vortex and the relaxation of magnetization are decoupled. When F^ decreases, the vortex is re- 
trapped by the background of remnant magnetization and they again form vortex polaron. This process repeats 
when FlU) increases in the opposite direction. Remarkably, after dissociation, decoupled vortices move in the 
periodic potential induced by magnetization which remains for some periods of time due to retardation after the 
decoupling. At this stage vortices oscillate with high frequencies determined by the Lorentz force at the moment 
of dissociation. We derive also the creep rate of vortices and show that magnetic moments suppress creep rate. 

PACS numbers: 74.25.Wx, 74.70.Dd, 74.25.Ha 
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I. INTRODUCTION 

To explain behavior of the critical current in Er borocarbide 
ErNiaBaC' ^ it was proposed in Ref. 6 that in magnetic su- 
perconductors vortex motion becomes suppressed due to their 
interaction with the nonuniform part of magnetic moment po- 
larization. Magnetic field in the vortex state is nonuniform and 
so is the polarization of magnetic moments. In the flux flow 
state, vortices, to move, need to change polarization as they 
are attracted to the polarization clouds induced at their previ- 
ous positions. Due to nonzero magnetic relaxation time vor- 
tex viscosity is enhanced and this enhancement increases with 
magnetic relaxation time. The vortex together with accom- 
panied nonuniform polarization form a vortex polaron at low 
vortex velocity. As bias current and Lorentz force increase so 
does vortex velocity and above some critical velocity vortices 
get rid of polarization and move at higher velocities defined 
by the usual Bardeen-Stephen drag coefficient. At this po- 
laron dissociation voltage jumps to a higher value and from 
the point of view of I-V characteristics the corresponding bias 
cuiTent may be identified as the critical cuiTent. When pinning 
due to defects is present, vortices stait to move with enhanced 
viscosity and low dissipation above the critical cuiTent, deter- 
mined by the pinning due to defects. Then above the dissoci- 
ation critical current they move with the velocity determined 
by the Bardeen-Stephen drag coefficient^. Thus effectively 
the critical current is the sum of standard pinning critical cur- 
rent and polaron dissociation critical cuiTent. The polaronic 
mechanism constraining vortex motion is effective not only in 
crystals, but in multilayered systems of superconductor and 
magnetic layers whose thickness is smaller or of the order of 
the superconducting penetration depth A; see Ref 7 

The polaronic mechanism is inherent to any magnetic su- 
perconductor, but the strength of magnetic coupling of vor- 
tex with magnetic moments depends crucially on the London 
penetration length which determines the nonuniform part of 
the magnetic field inside superconductors. The interaction 



drops as l/zt"* as /I increases and thus the polaronic mecha- 
nism is most effective in superconductors with low A's like 
borocarbides where A is of the order 500 A. Another impor- 
tant parameter is the magnetic susceptibility = dM/dH (M 
is the magnetization) at the operating applied magnetic field 
H. Finally the third crucial parameter in slowing down vortex 
motion is the relaxation time r of the magnetic moments. 

In this paper we consider vortex motion in the presence of 
the alternating in time Lorentz force and show that the pola- 
ronic effect results in suppression of vortex oscillations and, 
consequently, dissipation. We consider also thermally acti- 
vated vortex motion via the potential barrier inherent to vor- 
tex creep at low temperatures in the Andersen-Kim model. We 
show that creep rate is also suppressed by the polaronic effect 
resulting in lower dissipation at low temperatures. We predict 
that a strong ac Lorentz force with frequency oj ~ l/ixr) re- 
sults in vortex lattice oscillations with the frequency Q. '» co 
just after polaron dissociation due to retardation in response 
of magnetic moments to the inhomogeneous magnetic field of 
moving vortex lattice. 



II. GENERAL EQUATIONS 

We consider the system of vortices at the coordinates R,- - 
{xi,yi) interacting with the magnetic moments described by 
the magnetization M(r) in thin film with the thickness d much 
smaller than the London penetration length A. The Lagrangian 
X{Ri(0> M(r, t)} for the whole system reads 



£{Ri(t), M(r, 0) = XM{M(r, f)) + Xv{R,(f)) 
+Xi„,{M(r, f), R,(f)) + -CvARi) + Xf{J). 

The Lagrangian for the magnetic subsystem is 



£M{M{r,t)] = -d drM\r,t)l{2x) 



(1) 



(2) 



2 



The Lagrangian for the interaction between vortices and pin- 
ning potential due to quenched disorder reads 



X,{R,(0,r,) = -y f/(R, -r^), 



(3) 



where t/(R, -ry) is the pinning potential at r^. Further, -Cn lR;) 
is the vortex-vortex interaction and -£f{J) - E/ J ' ^i'^od/c is 
the Lagrangian due to Lorentz force in the presence of bias 
current density J, while Oo = hc/(2e) is the flux quantum. 
Here x is the magnetic susceptibility at the working external 
magnetic field. It describes response of magnetic moments to 
nonuniform component of the field induced by vortices. We 
describe the magnetic moments in the continuous approxima- 
tion via the magnetization M(r, t), because distance between 
spins is much smaller than the superconducting penetration 
length A. We consider the Ising spins with the moments along 
the vortex direction, as in Er borocarbide. We ignore the pair 
breaking effect of the magnetic moments. The interaction be- 
tween vortices at the coordinates R,-, and the magnetic mo- 
ments is determined by the term 



Xi„,{R;, M) = £/ drB(Ri - r, t)M{r, t). 



(4) 



Both the vortices and magnetization are described by a re- 
laxation dynamics characterized by the dissipation function 
:R{R,(f), M(r, 0) = + where 

KMim)] = \rd J drM\r), ^ndj] (5) 

Here r is the relaxation time for a single spin and rj is 
the Bardeen-Stephen drag coefficient per unit vortex length, 
rj - <l)y/(27r^^ c^p„) with p„ the normal resistivity slightly 
above Tc- Within the London approximation, the mag- 
netic field B(r) induced by vortices is given by B{r) - 
^/i2nf 2,- / t/kB(kexp[/k(r - R;)], k = (k^,ky) with 



B(k) 



On 



1 +^%2- 



(6) 



Note, that the superconducting penetration length here is the 
renormalized length given by the expression - Aj^(l -Ajix), 
where A^ describes magnetic field penetration in the absence 
of the magnetic moments The equation of motion for vor- 
tices is 



It gives 



d 6£ 6£ d-R „ 

-f — - TiT + — 
dt (5R,. 5Ki 5R,. 



m ^ ^i:vv{R,-,Rj} ^ gXint{R,-,M) 



(7) 



j 



5R; 



H-^/Fi + rfiCR,-,?), 



(8) 



where — OoJ/c is the Lorentz force. The last term is the 
Gaussian random force with the correlation function 

<rR,^(R,-, f)r«,v(R), ?')) = {2kBTr^d)6^6^,M - t'y (9) 



where T is temperature, kB is the Boltzmann constant and 
li,v - x,y. We assume that the dynamics of the magnetization 
is the dissipative one. 



dM{r, t) 
dt 



M{r, t) 



X 



-B{r) 



+ rM(r,f)- 



(10) 



From Eq. ( 10 1 we see that relaxation time of the magnetization 
measured experimentally in the crystal is xt- The last term is 
the Gaussian random force with the correlation function 

{TMir, tWMir', t')) = {2kBTT)6{r - r')5(t - f). (11) 



III. RESPONSE OF THE VORTEX LATTICE TO AN AC 
DRIVING CURRENT 

We consider the response of the vortex lattice to an ac driv- 
ing current in the flux flow state. In this case, the thermal fluc- 
tuations can be neglected. Meanwhile, the quenched disorder 
can also be neglected because the vortex motion quickly aver- 
ages out the disorder and the lattice ordering is improve d '^^^^ l 
In the absence of thermal fluctuations and quenched disorder, 
the vortex lattice moves as a whole and the interaction be- 
tween vortices cancels. The motion of the center of mass for 
the vortex lattice, u{t) , along the ;t:-axis is described by the 
equation [r - (x, y)] 



du 
dt 



d_ 
du 



drBix + u(t),y, t)M(r, t) 



+ Fr. 



(12) 



In the Fourier representation of the vortex lattice, we have 
coupled equations for the center of mass of the vortex lattice 
u(t) and components of the magnetization M(G, t): 



du V— I 



<Do exp[/G^-M] 



5M(G, f) 
' dt 



M(G,f) ^ Oo exp[-!Gj:M] 



X 



1 + A^G^ 



(13) 



(14) 



where n,, is the vortex density and G are the reciprocal vec- 
tors of vortex lattice. The coupling between magnetization 
and vortices drops as (1 -H G^A^y^ as G increases, thus it 
is sufficient to take the dominant lattice wave number Gi = 
(27r/fl,0), where a is the lattice constant^ ^. Here we consider 
a dense square lattice and n,, - l/a^, while GiA » 1. Renor- 
malizing time in unit of tx, length in unit of 1/Gi, force per 
unit vortex length in unit of rjlijGxx), we have equations for 
m(f) 



M{Gx,t)A^G\l{(^oX) and u{t): 

dim(t) - -[m{t) - exp[-;M(f)]], 

d,u = Fl -Im ^Fp exp(/M)m(f)] ■. 



(15) 



(16) 



with an ac Lorentz force F^ - F^^ sin(wf) and magnetic pin- 
ning force per unit vortex length Fp = (^X^^tHIjP'A^t]). Elim- 
inating m{t) we obtain equation for u{t): 

du r' , , , 

— - Fi - Fp dt sin[M(f) - u(t )] exp(f - t). (17) 
dt Jo 
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FIG. 1. (color online) (a) The dependence of effective viscosity rj^g, 
pinning strength ap, and (b) normalized dissipation power, D/Dg, on 
the driving frequency ai in the linear regime F^c (i>(i + Fp)- Here 
F„ = 20. 



was induced by the same lattice at previous positions and pre- 
vious moments of time. The important point for possible ap- 
plications is that in the presence of magnetic subsystem, the 
dissipation of the system is strongly reduced in linear regime 

For a stronger driven force amplitude, F^^ > Flc, in the 
hysteretic regime, we describe the system analytically in the 
adiabatic limit, w <«; 1 . At the time moment tc, when F^it) = 
Flc ~ 0.5Fp, polaron dissociation leaves magnetization and 
vortex lattice weakly coupled because lattice moves now with 
a high velocity. The magnetization component m(t) after that 
moment relaxes as m(t) - exp(-f + tc), and motion of vortex 
lattice is determined by the equation 



du 
dt 



Flc + Fp sin(M) exp(-f + tc). 



(21) 



The dc limit to -Q has been considered in Refs. |6]and|2l The 
behavior of the system depends on Fp. For Fp > 8 the sys- 
tem exhibits hysteretic behavior. Upon increasing the current 
from zero, at the critical velocity Vc ~ I the system jumps to 
the conventional Bardeen-Stephen branch by dissociation of 
polaron at the threshold force Fi - Fic « 0.5Fp. When one 
reduces the current in the Bardeen-Stephen branch, the sys- 
tem jumps to the branch with larger viscosity by formation of 
polaron at a retrapping current that is different from the disso- 
ciation current. 

In the ac current regime, for a low amplitude Fnc/[(y(l + 
Fp)] <sc 1, the vortex lattice oscillates, u - Re[Mflc exp(/wf)], 
with the amplitude 



Uac = Fac(;;7effW + Up) 



(18) 

Up = FpcJ-{o? + 1)-'. (19) 



For a low frequency w «c 1, the effect of magnetization is 
to renormalize the viscosity from 77 to tj^s - \ + Fp. In this 
polaron region, the magnetization follows vortex motion by 
formation of vortex polaron, as in the dc case w = 0, result- 
ing in enhancement of viscosity and suppression of ac dissi- 
pation. For a large frequency » 1, the effect of magne- 
tization is to introduce pinning potential Um - FpU^ j2 with 
strength Fp. In this case, the vortex lattice follows the driving 
force much faster than magnetization, which remains almost 
time independent. In the presence of magnetic moments the 
dissipation power of the whole system, averaged over time, 
D{a)) - {Fi(t)v(t)),, is reduced. In the linear region with vor- 
tex polaron, we obtain 



D 



(20) 



This dissipation power should be compared with that, 
DoF^c/2, without magnetic moments (at Fp - 0). For w «c 1 
we get DjDo = (1 + Fp)"' and for w » 1 we get D/Dq ~ 1. 
The frequency dependence of the effective viscosity, //eif , pin- 
ning strength Op and normalized dissipation power DjDo is 
shown in Fig. [T] 

The polarization of the magnetization results in periodic 
pinning potential with the periodicity of vortex lattice as it 



When t - tc > 1 the velocity of the vortex lattice oscillates 
with the frequency Q - Fic but oscillations relax on the time 
scale of unity. 



V ^ Fic+ Fp sm(Fict) exp(-f + tc). 



(22) 



These post dissociation oscillations are due to the motion of 
vortex lattice in periodic potential induced by remnant re- 
tarded magnetization. Results for arbitrary time t - tc aie 
shown in Fig|2] which are compared with the direct numer- 
ical solution of Eqs. ( 15 1 and (16 1. 

For arbitrary ll) we solve numerically Eqs. (15 1 and (16 1 



taking both the retardation and nonlinearity into account. We 
consider the interesting region Fp > 8, where the dissociation 
of vortex polaron is possible due to nonlinear effects at m > 1 . 
We take Fp - 20 in the following discussion. The hysteretic 
behavior of vortex lattice velocity vs. driving force is shown 
in Fig. [3] At frequencies a> < I, during the period of Fi(t), we 
see the following sequence of events: polaron formation near 
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FIG. 2. (color online) (a) v(t) and Fi(t) as a function of time obtained 
by numerical solution of Eqs. ijTSj and l[T6j with lo = 0.01, F^c = 20 
and Fp = 20.(b) Red linexlose up for the fast oscillations of velocity 
in (a). Blue line numerical solution of Eq. {21\ with Fic = 8.103. 
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FIG. 3. (color online) Dependence of the vortex velocity v(t) on the 
driving force Fi(t) = F.^^ sin(ajr) with = 20. Here Fp = 20. 



period of Fi{t) is significantly bigger than unity. This expres- 
sion for the frequency in original units reads Q » InFuc/aJ]- 
Such relation is anticipated for decoupled vortex moving in 
the pinning potential with periodicity a. In Fig. |5] we show 
normalized spectral function 5(Q) - \v^(Q.)\/FI^ for differ- 
ent driving force frequency at F-^^ - 20 [we denote v(Q) = 
J dtv(t)e\p(iD.t)]. The spectrum of oscillations consists of 
harmonics (2« + l)cj due to nonlinearity of the equation for u 
(n is integer) and amplitudes of these harmonics are enhanced 
near the frequency Q ^ Flc when oj < I. 

In the oscillation regime part of dissipation power, 
{v{t)Fi(t)), is transferred into oscillation power in the fre- 
quency region near Q - Fic. This part increases with the 
number of oscillations per half of the Lorentz force period. In 
Figs.|6]we show dependence of normalized dissipation power, 
D(F^c)/Do, as a function of F.^^. and the corresponding behav- 
ior of v(t) and F/,(f). For oj < I, dissipation increases step 
wise as F^^ increases. To understand the plateaus and steps in 
dissipation power, we plot v(f) and Fi(t) at several corre- 
sponding to plateaus in D/Dq. As shown in Fig.|6] the differ- 
ent plateaus correspond to different number of periods of the 
fast oscillations during the time n/oj. Hence, these plateaus 
are due to the commensurability in periods of fast oscillations 
and slow oscillations of the driving force. The sharp increase 
of dissipation power from one plateau to another is caused by 
changing of the period of fast oscillation. The behaviour in 
the region Fac/<^ » 1 is highly nonlinear, and numerical cal- 



low \Fi\ (interval of low vortex velocity), then polaron disso- 
ciation (velocity sharp increase) followed by region of vortex 
oscillations on the background of average high velocity, de- 
crease of velocity as the Lorentz force drops and vortex re- 
trapping (sharp drop in vortex velocity) and then dissociation 
again at negative -Flc (sharp drop in velocity). The results 
for behavior of vortex velocity in time, v(f), at F^c = 10 > Fj^ 
and different w are shown in Fig. |4] 

At all frequencies w < 1 we see post dissociation oscil- 
lations caused by motion of decoupled vortices with respect 
to periodic potential created by nonuniform magnetization in- 
duced by the same lattice just before decoupling (when ve- 
locity was still low) and frozen for some periods of time af- 
ter decoupling due to the retardation effect. This self-induced 
pinning, resulting from the retardation, and the amplitude of 
corresponding vortex oscillations reach maximum at oj ^ 1. 
In a rough approximation we describe them by the equation 



du 
dt 



Flc + FpMd sm(u - Ud), 



(23) 



assuming approximately constant m and F^ - Fi^ in the re- 
gions of maxima and minima of the Lorentz force. Here Ud 
and nid are the position of vortex lattice and the amplitude of 
magnetization at the moment of decoupling. This gives ap- 
proximate solution near extrema of Fi{t) at F.^^: 



v{t) ~ Flc + FpMd siniFict), 



(24) 



which provides rough estimate for the oscillation frequency, 
^ ~ Flc, when number of oscillations of velocity per the half 




K h h h h I 








FIG. 4. (color online) Time evolution of the vortex velocity, v(t), and 
of magnetization \m{t)\ in the presence of ac driving force Fi(t) = 
sin(<i;/) at several frequencies ai. (a): cj = 0.1, (b): a) = 1, (c): 
OJ = 10. We take F,, = 20 and F„ = 20. 
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FIG. 5. (color online) Normalized spectral function S(C1) of vortex 
lattice oscillations at different frequencies of the driving force a> = 
0.1,0.25, 1, 10 at the amplitude F^c = 20 and F,, = 20. 



culations become important tool to understand the dissipation 
in the system. 

For high frequency w » 1 it is impossible to accommo- 
date fast oscillation with frequency F^^ during the period of 
the external driving frequency, and step structure disappears, 
as shown in Fig. [7] For high frequency co the magnetization 
m(t) is approximately periodic in space (as induced by peri- 
odic vortex lattice) and oscillates in time with the amplitude 
which decreases with ll) as shown in Fig.|4] 

Hence, at low driving force amplitudes, F-^^ < Q.5Fj,, the 
polaronic effect renormalizes the viscosity at low oj and in- 
troduces pinning at high co. At a long magnetic relaxation 
time, ojTx » 1, the elastic constant of the magnetic pin- 
ning force per unit vortex length in dimension form becomes 
ap = FpTiKGiXT) = <l>,^/(27r2.i'*Gi), while at small utx « 1 
the vortex viscosity is enhanced from rf to 77(1 +Fp). At higher 
values of reduction of dissipation is lower. In Fig. |7]we 
show the normalized dissipation power as a function of F.^^ 
and Lo. Due to commensurability effects as shown in Fig. |6] 
Z)o/Z) oscillates with (±> in the region (±> \ when F > F 
as shown in Fig. |7](b). For a large u or Fac < Flc, Dq/D 
increases smoothly with co. 

To understand significance of polaron pinning we com- 
pare the elastic constant of pinning potential for vortex o-p - 
FpTjl (jx) with the maximum possible Op for a single vortex 
pinning. Such a maximum, Uc ~ H^/(8n) - ^^^/i64n^A^^^) 
is reached in the case of columnar defects produced by irradi- 
ation of crystal by heavy ions''*, where we have assumed that 



the size of defects is ^, while He the thermodynamic critical 
field. 

For ErNi2B2C at temperature 3 K, A = 500 A, ^ = 135 A, 
X ~ 0.05, p„ = 10"^ Q ■ m and we estimate ap/a^ x 0.4.E3[ni 
Thus magnetic pinning is not drastically lower than maximum 
possible one. 

To observe hysteretic behavior and high frequency oscil- 
lations of voltage following the dissociation, one needs to ex- 
ceed the critical current Jc corresponding to Flc. In ErNi2B2C 
this critical current at w — > is quite high, about 10^ A/cm^ 
at low temperatures for the field 0. 1 T along the a or b crystal 
axis. It increases with oj as shown in Fig. [3] However, Jc pro- 
portional to sin^ a becomes smaller , as field approaches the 
c-axis, due to renormalization of x x sin^ a. Here a is the 
angle between the applied magnetic field and the c axis. 

The important parameter which determines the dynamics of 
polaron vortices in magnetic superconductors is 

2^2 2 

Fn-^^- (25) 

It depends crucially on the magnetization relaxation time x^, 
measured experimentally, and the magnetic susceptibility x 
which are only known for the most studied Er borocarbide. 
In this crystal according to Mossbauer measurement^ ;^fT ^ 
5 X 10"'" s at r = 5 K and we anticipate longer time at lower 
temperatures. For x - ^ ^ 10"^ and p„ - 5 /^Q-cm we get 
Fp > 120. In Tm borocarbide TmNi2B2C we anticipate big- 
ger at temperatures above Neel temperature about 1 K due to 
larger density of almost free spins, but t for Tm borocarbide 
is unknown. 

As a consequence of polaron dissociation at driving force 
with high amplitude and not very high frequency vortex lattice 
velocity oscillates and generates ac voltage with frequency of 
the driving force and also high frequency voltage whose fre- 
quency depends on the driving force amplitude. These high 
frequencies are about 10-20 times bigger than driving force 
frequency and conversion of the power to these frequencies 
is quite efficient at w < 1 and F-^^ ~ 20. This implies that 
the vortex motion in magnetic superconductors may be used 
to generate high frequency electromagnetic oscillations whose 
frequency depends on the amplitude of the driving current. 



IV. VORTEX CREEP 

Next we discuss how polarization of the magnetic moments 
affects the vortex creep. We consider motion of a single vortex 
from one pinning center at ro to another pinning center at ri 
due to the thermal fluctuations, as schematically shown in the 
inset ofFig.|8] We assume that Uiro) > U{ri). The behavior 
of vortex is characterized by a rate for thermal escape of vor- 
tex from a metastable potential. In the absence of the magne- 
tization, the rate can be calculate using the standard Kramers' 
theory'**. The magnetization with relaxation dynamics intro- 
duce memory effect for the vortex crossing the barrier Here, 
we investigate the effect of retardation on the escape rate. 

Without a Lorentz force, the equation of motion for vortex 
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FIG. 6. (color online) The normalized dissipation £)(Fac) as a function of the amplitude of the driving force and time evolution of v(t) and 
Flit), at several values of F^c, corresponding to the different plateaus in the dissipation power D/Dg in (a). 



kxMik) Oo exp(il,r) 
^2k2 + 1 



+ r«, (26) 



r 

',r = -drU + dk- 
r5,M(^) = - (— - ^ ^^5,^ j + Tm. (27) 



Main contribution to integral in Eq. ( 26 1 comes from the re- 
gion k < I /A. We assume that the distance between pinning 
potentials is much smaller than A, i.e. |ri - ro| <s: A. Then we 



can expand in k^r in Eqs. (26 1 and (27 1. We then obtain the 



following equation of motion for vortex 

d,r = -drU + a(Mi -r) + T,, 
XTdiMi = - (Ml -r) + ri, 

with the correlation functions 

{Y,(t)TAt'))^2T6(t-t'), 
{r,(t)ri(t'))^2TT6(t-t')/a, 



(28) 
(29) 

(30) 
(31) 



where we have introduced the following dimensionless units: 
r and Mi is in unit of ^, U is in unit of Uq, time is in unit 



of T]^^ IU{), temperature is in unit of U[)d/kB- The strength of 
coupling between magnetization and vortex is a = A^^/Uo- 
Here 



Op 
2nA 



r ik,M(k) 
J'^'a^V 



J (2;r)2 ^2k2 + 1 4^^4 [ \ ^ j 



(32) 



■(33) 



From Eq. (27), we obtain for Mi 



Jo 



dAt') ■ 



ri(f') 



dt', (34) 



Mi(f) = r(f)-e- 

with r' = ;ifT. Substituting Eq. ( |34| ) into Eq. ( |28] l, we obtain 

r /:(f - t')d,r{t')dt' = + r, (35) 

Jo 



with 



K{t - t') - aexp[-|f - f'l/r'] + 5{t - t'). 
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FIG. 7. (color online) Dependence of the normalized dissipation 
power D/Dq on F^c in (a) and on a> in (b). Here Fp = 20. 
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FIG. 8. (color online) Renormalized rate R{t')/R{t' = 0) as function 
of t' at several a. Here we take a>ij = 1 . Inset is a schematic view of 
thermally activated vortex escaping from a pinning potential. 



Jo r' 



Note that the nonlocal dissipation and noise still satisfy 
fluctuation-dissipation theorem 



(r'(f)r'(f')) = 2TK{t - 1'). 



Escaping rate for viscosity with retardation eff'ect was calcu- 
lated by Hanggi and Mojitabai in Ref.[19] We use their results 
for our case of the viscosity containing parts with and without 



retardation. We obtain the expression 

^(•^') = exp[-Af//r], 



(36) 

^ (37) 



where wo - ^/^jU(rQ), cjh - y/-d'^U(rb), and AU - U(rh) - 
U(ro). Here is the position of the barrier top. The results 
for R(t')/R{t' = 0) are shown in Fig. |8] The rate decreases 
with r'. Equation (36i is valid if R{t')t' <«; 1. In the limit 



of strong interaction of voitex with magnetic system, a > uj^, 
this expression is valid at any r because ji drops with t as 
1/t^. In this case the barrier remains the same as in the system 
without magnetic moments (at a - 0), but pre-exponential 
factor diminishes as t increases due to interaction of vortex 
with magnetic moments. 

For weak interaction, a < co-^^, the rate has nonzero limit at 
a large r, and the condition 7?(t')t' <k 1 may be violated at a 
large r. When it occurs, the barrier is renormalized from AU 
to AU + a{rh - ro)^/2 as magnetization cannot follow escaping 
vortex. Pre-exponential remains the same as the bare vortex 
system with a - 0. 

In both cases the rate decreases due to the polaronic eff'ect. 
Eff'ect depends strongly on the value a which is unknown for 
ErNi2B2C at this moment because the characteristics of the 
pinning centers are unknown. Measurements of the relaxation 
of magnetization lilnM/iilnf above 6 K, where the polaron 
mechanism is absent, may reveal the pinning barriers in the 
crystals. Then one can predict suppression of the creep rate at 
lower temperatures. 



V. CONCLUSION 

In conclusion, the polaronic effect in magnetic supercon- 
ductors (crystals and multilayered structures) suppresses dis- 
sipation for ac currents. This eff'ect depends strongly on the 
parameter wry. When this parameter is large, magnetic mo- 
ments result in additional strong pinning due to the polaronic 
effect. Creep rate at low temperatures is reduced as well. Thus 
magnetic superconductors may be useful for applications, es- 
pecially for devices with ac cuiTents. We predict the gener- 
ation of voltage oscillations with high frequency il » w in 
the presence of strong ac driving Lorentz force with moderate 
frequency cj ~ 1/(tx)- This effect is caused by pinning of 
vortex lattice on retarded periodic magnetization induced by 
the same vortex lattice in previous moments of time. 
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